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1. Introduction 

Improved staggered fermion formulations are widely used to study QCD at non-zero temper- 
atures and densities, see e.g. Ref. ||] for recent reviews. The main reason for this is the fact 
that they preserve a part of the chiral symmetry of the continuum QCD. This makes the numerical 
simulations relatively inexpensive because due to absence of an additive mass renormalization the 
Dirac operator is bounded from below, and also allows one to study the chiral aspects of the finite 
temperature transition. However, there are at least two problems with staggered fermion formu- 
lation. The first one is the validity of the rooting procedure, i.e. the way to reduce the number 
of tastes from four to one and the other is breaking of the taste symmetry at finite lattice spacing. 
For the discussion of the validity of rooted staggered fermions see Refs. [^, |]]. To reduce the 
taste violations smeared links are used in the staggered Dirac operator and different staggered for- 
mulations, like p4, asqtad and stout differ in the choice of the smeared gauge links. The smeared 
gauge links in the p4 and asqtad actions are linear combinations of single links and different staples 
[]|, ^] and therefore are not elements of SU(3) group. It is known that projecting the smeared gauge 
fields onto the SU(3) group greatly improves the taste symmetry [f7|]. The stout [^] action and the 
HISQ action considered here implement the projection of the gauge field onto SU(3) (or simply 
U(3)) group and thus achieve better taste symmetry at a given lattice spacing. For studying QCD at 
high temperature it is important to use discretization schemes which improve the quark dispersion 
relation, thus eliminating the tree level 0(a 2 ) lattice artifacts in thermodynamic quantities. The p4 
and asqtad actions implement this improvement by introducing 3-link terms in the staggered Dirac 
operator. In this contribution we report on exploratory studies of QCD thermodynamics with the 
HISQ action which removes tree level €?(a 2 ) lattice artifacts as well as has projected smeared links 
that greatly improve the taste symmetry. We also compare our results with previous ones obtained 
with the asqtad, p4 and stout actions ^ [i0|, 12, 13 J. 



2. Action and run parameters 

The Highly Improved Staggered Quark (HISQ) action developed by the HPQCD/UKQCD 



collaboration [14] reduces taste symmetry breaking and decreases the splitting between different 
pion tastes by a factor of about three compared to the asqtad action. The net result, as recent studies 
show [jig], is that a HISQ ensemble at lattice spacing a has scaling violations comparable to ones 
in an asqtad ensemble at lattice spacing 2/3a. 

In this exploratory study we used the HISQ action in the fermion sector and tree-level Symanzik 
improved gauge action without the tadpole improvement. The strange quark mass m s was set to its 



physical value setting the quantity y 2m 2 K — m 2 l = m^- ~ ^2Bm s to the physical value 686.57 MeV. 
Two sets of ensembles have been generated along the two lines of constant physics (LCP): m/ = 
0.2ra s and m/ = 0.05ra. s . In the first set runs were performed on 16 3 x 32 lattices at zero temperature 
and 16 3 x 6 at finite temperature, and in the second set on 32 4 and 32 3 x 8 lattices, correspondingly. 
The parameters and statistics of these runs are summarized in Table ||. The molecular dynamics 
(MD) trajectories have length of 1 time unit (TU) and the measurements were performed every 5 
TUs at zero and 10 TUs at finite temperature. Normally, 300 TUs were discarded for equilibration. 
The lattice spacing has been determined by measuring the static quark-anti-quark potential and us- 



2 



First results on QCD thermodynamics with HISQ action 



Alexei Bazavov 





N t = 6, 
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Table 1: The parameters of the numerical simulations: gauge coupling, strange quark mass and the number 
of time units (TU), i.e. the number of MD trajectories for each run. Here TU, stands for the number of 
time units for zero-temperature runs, while TU, T is the number of time units for finite-temperature runs. 
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Table 2: Ranges of masses (in MeV) of several hadrons for the two sets of ensembles. 

ing the Sommer scale ro = 0.469 fm. The masses of several hadrons measured on zero-temperature 
ensembles fall into ranges summarized in Table |[ 

As mentioned above the lattice artifacts for the HISQ action are significantly reduced com- 
pared to the asqtad action. The taste violations are strongest in the pseudo-scalar meson sector. 
Therefore we studied the splitting of pseudo-scalar meson masses in the 8 different multiplets. The 
results are shown in Fig. [T] and compared to the stout results. The splittings are 2 to 3 times smaller 
than in the calculations with the asqtad action and also somewhat smaller than for the stout ac- 
tion. The smaller taste violations also result in better scaling of other hadron masses. In Fig. [j] we 
also show the kaon, rho-meson and nucleon masses. As one can see for the small, almost physi- 
cal, value of the light quark mass we get reasonable agreement between the lattice results and the 
experimental values. 

3. Deconfinement transition 



The deconfinement transition is usually studied using the renormalized Polyakov loop [12. 
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Figure 1: The splitting between different pion multiplets (left) and hadron masses (right) calculated for 
the HISQ action. In the right panel squares indicate ensembles along the 0.2m s LCP, and circles along the 
0.05m v LCP. 



17], pTSp . It is related to the free energy of a static quark anti-quark pair at infinite separation Foo(r) 

L ren (T) = exp(-Feo(r)/(2r)), (3.1) 



obtained from the bare Polyakov loop as 



Lren{T)=z{P)^L bare {p)= Z {P) 



1 ATt-1 * 

T Tr EI ^o(*0,*) 
3 *o=0 , 



(3.2) 



Here the multiplicative renormalization constant z(j3) is related to the additive normalization of the 
potential c(j8) as z(/3) = exp(— c(j8)/2). To make the comparison with the stout results easy here 
we chose the normalization constant c(j3) such that the potential is zero at distance r = r^. This 
choice is different from the one used in Refs. 0, |TJ]. Our results for the quark mass m q = 0.05m s 
on 32 3 x 8 lattices are shown in Fig. || and compared to the stout and p4 calculations. As one 
can see the HISQ calculations agree reasonably well with the stout results if the scale is set by 
ro in the stout calculations 1 . On the other hand, the renormalized Polyakov loop calculated with 
the p4 action is considerably smaller at low temperatures. At temperatures T > 200 MeV we see 
good agreement for different actions. The fluctuation of strangeness is also a good indicator for the 
deconnnement transition. It can be defined as the second derivative of the free energy density with 
respect to the strange quark chemical potential 



Xs(T) 1 d 2 \nZ{T^ s ) 



T2 T^V d(ji s /T) 



(3.3) 



At low temperatures strangeness is carried by massive hadrons and therefore strangeness fluctua- 
tions are suppressed. At high temperatures strangeness is carried by quarks and the effect of the 
non-zero strange quark mass is small. Therefore, in the transition region the strangeness fluctuation 
rises and eventually reaches a value close to that of an ideal quark gas. Our numerical results for 



'in what follows we will show the stout results using the temperature scale set by ro instead of the kaon decay 
constant f%. We used the published values of and f% in Refs. | p^ , |TJ] to convert the two scales. 
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Figure 2: The renormalized Polyakov loop as function of the temperature calculated for the HISQ, stout and 
p4 actions at the physical value of light quark masses. 

strangeness fluctuations are shown in Fig. || for two values of the light quark masses, m q = 0.05m v 
and Q.2m s . In the case of the heavier quark mass we compare our results with the ones obtained 
with the asqtad action for m q = 0.2m s . As one can see the strangeness fluctuations are larger for 
the HISQ action at low temperatures, T < 210 MeV. In other words, the transition region in the 
HISQ calculation shifts toward smaller temperatures. This behavior is in fact expected. Due to 
smaller taste symmetry violation pseudo-scalar masses as well as other hadron masses are smaller 
and therefore strangeness fluctuations are larger. For the smaller (physical) quark mass we compare 
our calculations with the results obtained using the p4 action [[K], [n]] and the stout action [12, |l3| ]. 



At low temperatures, T < 200 MeV the HISQ results are significantly larger than the p4 results but 
are in good agreement with the stout results. At high temperatures, T > 200 MeV the strangeness 
fluctuations calculated with the HISQ action are in reasonable agreement with the p4 results as well 
as the N t = 12 stout results. We also see that in this temperature region the stout results show some 
cutoff (N t ) dependence. This is due to the fact that the tree level 0{a 2 ) lattice artifacts are not 
removed in the stout action. In the low temperature region we expect that the Hadron Resonance 
Gas (HRG) model gives a reasonably good description of the thermodynamic quantities, including 
strangeness fluctuations. Therefore in Fig. || we show the prediction of the HRG model. As one 
can see the lattice results fall below the HRG model results. The discrepancy between the HRG 
model and the lattice data becomes larger at smaller temperatures, although the model is expected 
to be more reliable at smaller temperatures. This is presumably due to taste violation, especially in 
the pseudo-scalar meson sector. Due to non-negligible splitting in the pseudo-scalar meson mass 
the contribution of kaons to strangeness fluctuations appears to be smaller than expected in the 
continuum. 



4. The chiral transition 

In the limit of zero light quark masses QCD has a chiral symmetry and the finite temperature 
transition is a true phase transition. The order parameter for this transition is the light chiral con- 
densate (Y¥)i- However, even at finite values of the quark mass the chiral condensate will show 
a rapid change in the transition region indicating an effective restoration of the chiral symmetry. 
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Figure 3: The strangeness fluctuations % S (T) calculated for m q = §.§5m s (left) and m q — Q.2m s (right). For 
the physical quark mass m q = 0.05m s we compare our results with previous calculations performed with the 
stout and p4 actions. The solid line in the left plot shows the prediction of the HRG model. 



Since the chiral condensate has an additive ultraviolet renormalization we consider the subtracted 



chiral condensate [17, 18] 

(yy)/.r-|(wKr 
^ (r) -(wko-^(wko- (41) 

Here the subscripts / and s refer to light and strange quark condensates respectively, while the 
subscripts and % refer to the zero and finite temperature cases. In Fig. || we show the renormalized 
chiral condensate calculated with the HISQ action and compare it with results obtained with the 
stout action [ |13| ] as well as with the p4 results [10, [Tl| ]. Our results agree resonably well with 
the N t = 8 stout results and the agreement between the HISQ and stout results is even better if one 
considers N t = 10 and = 12. On the other hand, the subtracted chiral condensate is considerably 
smaller than for the p4 action. This is presumably due to the larger taste violating effects in the p4 
case. We also considered the fluctuations of the chiral condensate, which is just the disconnected 
part of the chiral susceptibility Xdisc- The corresponding results are shown in Fig. ^. For true phase 
transition the chiral susceptibility should show a peak at the critical temperature. The numerical 
results indeed show some peak-like structure. However, as discussed in Ref. [ JT9| ] it is not easy 
to differentiate between the increase in the chiral condensate due to nearby critical point and the 
effect of the Goldstone modes. Certainly, much more detailed studies are needed here. 



5. Conclusions 



In this contribution we discussed first results on QCD thermodynamics with the HISQ action. 
The use of this action allows one to minimize cutoff effects in thermodynamic quantities due to 
lattice artifacts in the quark dispersion relation and taste symmerty breaking. In fact, the HISQ 
action gives the smallest taste violation in the pseudo-scalar meson sector if measured though the 
mass splittings. We have calculated the renormalized Polyakov loop, the chiral condensate and 
the strangeness fluctuations on N t = 6 and N t = 8 lattices. We find good agreement between our 
results and the ones obtained with the stout action in the low temperature region, T < 200 MeV. At 
the same time present results disagree with the p4 results. In the high temperature region, where 
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Figure 4: The subtracted chiral condensate (left) and the disconnected part of the chiral susceptibility (right). 
For the subtracted chiral condensate we also compare our results with calculations performed for stout and 
p4 actions. Note that in the right panel we show the disconnected chiral susceptibility for the HISQ action 
for two LCPs: Q.2m s and 0.05ra s . 

there are visible cutoff effects in the stout calculations the HISQ results agree resonably well with 
the p4 results as well as with the stout results obtained on = 12. 

References 



II 

[2 
13 
14 
15 
16 
17 
18 
19 
[10 

[11 
[12 

[13 
[14; 
[15 
[16 
[17 
[18 
[19 



C. E. DeTar, PoS LATTICE2008, 001 (2008) 
P. Petreczky, arXiv:0908.1917 [hep-ph]. 
S. R. Sharpe, PoS LAT2006, 022 (2006) 
M. Creutz, PoS LAT2007, 007 (2007) 

F. Karsch, E. Laermann and A. Peikert, Nucl. Phys. B 605, 579 (2001) 

K. Orginos, D. Toussaint and R. L. Sugar [MILC Collaboration], Phys. Rev. D 60, 054503 (1999) 

A. Hasenfratz, Nucl. Phys. Proc. Suppl. 119, 131 (2003) 

Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, JHEP 0601, 089 (2006) 

C. Bernard et al. [MILC Collaboration], Phys. Rev. D 71, 034504 (2005) 

M. Cheng et. al. (RBC-Bielefeld Coll.), arXiV:091 1.2215 

P. Petreczky, these proceedings 

Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, Phys. Lett. B 643, 46 (2006) 

Y. Aoki, S. Borsanyi, S. Durr, Z. Fodor, S. D. Katz, S. Krieg and K. K. Szabo, JHEP 0906, 088 (2009) 

E. Follana et al. [HPQCD collaboration and UKQCD collaboration], Phys. Rev. D 75 (2007) 054502 

A. Bazavov et al. [MILC Collaboration], PoS LAT2009, 123 (2009) 

O. Kaczmarek, F. Karsch, P. Petreczky and F. Zantow, Phys. Lett. B 543, 41 (2002) 

M. Cheng et al., Phys. Rev. D 77, 01451 1 (2008) 

A. Bazavov et al., Phys. Rev. D 80, 014504 (2009) 

S. Ejiri et al, arXiv:0909.5122 [hep-lat]. 



7 



